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ABSTRACT 


Аг implicit numerical aodel for two-áinsnsional hydrody- 
nanic flow in coastal seas by Leendertse (1967), as modified 
by Hart (1976), was applied to Montsrey Bay. The model was 
tested against available water-levei and currant observa- 
tions. The responses of Monterey Bay to tidal forcing and 
steady-state winds were simulated. Under tidal forcing it 
was found to provide reasonable estimates of sea-surface 
elevations. Currants were not well predicted, indicating 
that other mechanisms such as wind, density stratification, 
and oceanic currents generally dominate the forcing of the 
circulation in Monterey Bay. The model in its present form 
was found to be potentially suitabls for providing real-time 
tide correctors during a hydrographic survey, achieving an 


RMS error of 4.5 cm in predicting s22-surfac2e elevations. 
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А. PURPOSE 


This study grew out of a desire to extend *idal data 
observed at a few iocations to the sentire area of a hydro- 
graphic survey. With sea-surface elevations modeled over 
the whole field, survey depths may be corrected automati- 
cally by subtracting realistic values of the surface's 
variaticn from datum at any point, at any time. 

Tidal zoning to obtain sea-surface elevations is, at 
present, a subjective affair requiring numerous water-level 
stations to indicate the progress of a tidal wave into an 
inlet. The pattern of propagation at points distant from 
the observing stations is inferred only qualitatively from 
bathymetry. Correctors are determined by defining zones 
graphically and computing appropriate phase and amplitude 
adjustments f3r each zone to apply to tidal values observed 
at the reference stations. This process dces not provide a 
continuum of correctors. It requires subjective judgment 
and considerable experience to échiave adequate results. 
The analysis is typically performed weli arter the survey, 
when it is too late to use depth data corrected for observed 
tides +o rrovide cross-checks on tha positional accuracy of 
the data. Errors that might have baen detected and 
corrected during the survey may pass unnoticed until an 
expensive return to the survey area or downward classifica- 
tion of the survey is necessary. 

Use in the field of a two-dimensional, numerical model 
for circulation and sea-surface elzvation wouid alleviate 
these difficulties. Such a model must be simple and flex- 


ible if it is to be applied in real or near-real «ime on 


10 





ield-type microprocesso-zs. In coastal seas and inlets, it 
should be able to provide sufficiently accurate surface 
elevations. 

To test this concept, the two-dimensional, hydrodynamic 
model of Leendertse (1967), as modified by Hart (1976) and 
during this study, was applied to Monterey Bay. The nodel 
had previously b¢én used with some success in shallow 
coastal seas and estuaries, but not in an area with such a 


SHEMMONTSLDSY canyon. 


л 


dominant bathymetric feature a 


ct 
Q 


Although no attempt was made compar2 this model to any of 
the broad spectrum of models in use throughout the oceano- 
graphic and coastal-engineering communities [ Tracor, Inc., 
1971], its relative simplicity, ease of implementation, 

NE UDlIity, and accurate output are all important to 1+5 
potential usefulness as a tidal zoning system during field 
operations. 

Momaddtcional purpose Sf this study is to incorporate 
large-scale non--:idal fozces into the model to =xplor2 their 
effects on the circulation and sea-surface elevation of а 
coastal body of water. Various investigators have 
Suggested, in fact, that tidél forces are overridden in 
their efrect on the circulation of Monterey Bay by the 
influence of »sffshore currents and atmospheric conditions. 
The potential significance of such factors, and the ability 
tO incorporates roal--zme observations of them into the 
model, are also importan* to the model's usefulness és a 


па) zoning system. 


B. HISTORY OF THE MODEL 


The original versicn of the numerical model used durir 
“his study was described by Leendertse (1967). His 
"“multioperational"” finite-difference scheme, usina both 


Ира Сз and axplicit techniques to s 


0 


lve the equations of 
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Би Шиљек оп, provided advantages la computaticnal stability 
and efficiency over the explicit moials then current [{ Hart, 
1976]. In particular, the model remains stable regaráless 
of the time step used; in a relatively deep embayment such 
as Monterey Bay the investigator is not restricted to time 
steps of the order of seconds as is common with explicit 
models such as that used by Lazanoff (1971) in his study of 
<һе bay. 

The model has been widely applied, both in small harbors 
[Leendertse, 1967; Leendertse and Lu, 1975; Chiang and Lee, 
19821 and in coastal seas [Leendertse, 1967; Har+, 1976; 
Spaulding and Beauchamp, 1983]. In past applications, 
length scales ranged to 290 km and iepths ranged to 100 m. 
This study extends the model +o a much desper area with 


pronounced vertical relief. 


С. CIRCULATION STUDIES OF MONTEREY BAY 


Monterey Bay is a relatively large (16 by 42 km), nearly 
symmetric embayment on the central coast of Califorria. Its 
most notatle bathymetric feature is the Monterey canyon, 
which curves into the bay from the southwest, severing the 
continental shelf. Within the bay proper, depths rise from 
750 m at the seaward end of the canyon to an average 55 m on 
the shelf. 

The bay presents to tha Pacific Ocean one uninterrupted 
open boundary, some 36 km long. Consequently, oceanic tides 
and currents and offshore atmospheric 2ffects are primary 
па forces of circulation within the bay itself.  Locel 
winds and seasonal river runoff may have some effects, espe- 
cially in the shallower portions of the bay to the north and 
scuth. The relative importance of these various influences 


is not well-understocd. 


12 





"орто тәрони Битиоц5 ”еәтуү Ápnas 9u3 јо deg LL sındra 
















































































| nn ი” — e o "” ...--------------___ე_.___--__ .«იი.”.ი–/წლშწეწ/წ'””””–  '–-._ _.__-ეებეაებეეაეაენეაბნბსასაეა..ი”,””-–ჟ9ოიც#–= –"!»”"""ო"ზ""შ'"ო"ჟ_.___=_..–_ ee ee _ღ_._._-__.__ უელ «=:.=- шш с — 
„с! SU 
ERIT EE Lc аы. Ийт: Te TU SE Eum "mfi Ті Aet I er u. FU mb I: EI H . > 90 
t 
+. п | 
уз % 9 შირ! Š deca ; вө}, с. аре иону; МІ SONIANNOS | 
“Фо, сан ^, წმ ҮП An С 0027 | || 
5 H E MISA SERI : 
| Е та зада озеђипр “ aa — E. єс»! | 
|Е = 6691 ү orez | 
Q . E 7 
: 7 990? } 
Anc EP > [474] - f ' 
| i : а + za ' f 6412 care 062 | 
a ‘Ore = 
j | qma m 4 2 = 2081 m = 
1 па E А dic CMS 
H ETT (er, — mt | { 
ЕБЕ N 
|| DE. NE ove 1612 \ 
в 
ti мл exi 22 oes! «0! თ DAN A x 
5 1981 -. инг — ғ \ 
$ yi ს ზი 5 
| | Тар j Ju. 
Š ; % (44.1) i —L- = — 
B ' ga ама | — ог | Z - I = 77 002 і 
-- 
_ i — 5191 bea ' | | тр — an ale a BF 
H y 2 i 2% ^ 1601 4 “o, = ee 4 ДЕ c ლი 
H E Ets t | } g zz è 
Е etu ` ТАР t s8. — | t 
Је ' აი s Ir а s ed ғо 
СН > МО вор eca < T tpe nmi... %L. № cove sr HPER 
a : с". | — I o, r | | j £ 
Op: . 25 \ шағала 
| туе) то Pa ~ j 2 ШЫ” ) 1%! 1 ЕР 
' шек 7760: ' 65/! „ი e 
| ао aS а т Y ‚ә ' %0% ( 
^ ey * | v ი)? en — e 1 
| LE № кыа +” Be G | 
e бк 96 а-а” ы i 0612 
sof usg | i] ve ў LT — о s d სა.“ К > 5 
: 4 ле war"! სი 975 ს ას თ шш C 
| H о 7261 1 > А T იის, - 
5 5 საა 99 ы. 1 
4 te, ^^. n8! -+05 J po қ | оса! д 
H ^ С 4 Set коз *-H 
è ღლ» eR 99 А. "i $ A ) f | + | 
% `` M M 09; - p ~ I ' 
а x У t 7 сз ve Е і 
3 დ. 4 E 5 бе v ! Б г у e 
I 4 ptt ч о УХ Nr One H 
| Н 5 en N m лон esa L 
1 თ у t€ Tex перељем ` \ 069! H 
| L ah - B 9C 79 `& EMN : | | 
= y) (789 * t 
' | . 30-71 = 4%, = em $ + 
E > E 1 ev |е / ` Н 
| E гиен > Е ce აა. ‚DIN ; 
| BE —— E А E Шы - — — 3 2 
с 5 4 — — — ~e A ⸗ |] > 
| E fu EN af pu LM — ee ec , | 454 Е 
` № - tt: : ta 8 z 
| ( სათესი рге Мо, ја — 54 60“ Ts : eoe _პ (0009 : 
д = Se (წ % ; . = 
კ t Nuez “\ lo, (9 : M 1 % сті 0с61 + г” © 
1 PU NA, + C P TL %%02494! £06€ 
= tv ` 1.2. <‘ | a aa 27? E 
ns 07 бе о ^" \ ან ! ^. . 8002 een? 54 
а ола , s რა ` tu. А = 
р — 0o ი. 004% sce ` 3 
$99 у weis nisus Oe > em DAT & ~, “о 
ә ғын vn жн та E 920 x 
o 222 რია 
| Hr mar japon p pmpn peg перо E түп це шт pnm SETE ТІ аг 
eccl ЕСІ Fl 


— — —À —ñ— — — — — — — — —— — — — — —s — — ——— — —— — — — —— ee — —— —— — — —— — ee ce ee ee аа 


13 





—————————————M—Ó—— it Жош. "Ры ы ““=”“”რწნ“”“=ო" —— э 






| 
I а = ) 
| SI a 75 RN RR — | 
5495 X ER 
BEEN DAN AN w | 
| SEER INVE a AO N | 
RENNER —— 
| ХХХ, QS УМУ M NN УШ Ж OR 
RUE GT ER — 
RY RR ს. ბა. LN 9 — EL X 
ba 2 5% N NN TAY 242 ER 
2% UY — * | 
N IR 
| “4 NE: А ER Am AN > | 
1 
г. AUS IT de — 
| 3 RR BR იას A Hi n 2 77 | 
| 4 У УУ ER — RER Y PLU HN EEE C | 
| ს eqs OR w A RN У MEER Ит | 
| NS = y VI ЖОРУ H | 
бо N | 8 УШ EHRT | 
| 1. ER | 
RE 224 | 
| ү ^0 ( у, > | 
Km =" | 
| ) 
| 20 | 
| 
0 | 
J 





Figure 1.2 Bathymetry of Monterey Bay Viewed from Southwest. 


Three oceanographic "seasons" for this portion of the 
coast were first described by Skogsberg (1936). From 
November through February the Davidson Current flows north- 
ward along the coast in conjunction with southerly or weak 
northward winds and the cnshore transport of surface water 
From March through August a period of upwelling is accompa- 
nied by the southward flowing California Curre2nt, strong 
northwest winds, and offshore transport of surface water. 

In September and October the oceanic period is marked by 
relative calm and an increase in surface temperatures. 

All available data concerning circulation within the bay 
proper were summarized and analyzed in 1973 as part of a 
Bor oceanographic study [ Scott, 1973]. Normal circulation 
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Ше “9108 290956 northward through thea bay, with a smail gyre 
erning in the Southern bight. A more recent analysis by 
Broenkow and Smethie (1978) also concluded that flow is 
generally northward, with water entering primarily along “he 
axis cf the canyon and having a rasidence tine of 2 to 14 
days during upwelling periods. They also suggest chat a 
volume of 109 m3 pumped into and out of the bay by internal 
tidal mixing may be responsible for the frequent pres O 
cool, nutrient-rich waters at tne head of the Montere 
Canyon near Moss Landing. 

Most circulation studies of th2 bay proper have relied 
primarily on temperature and salinity data collected at 
oceanographic stations. Although drift cards and similar 
devices have been used to map surface currents, long-term 
current-meter observations have not been available until the 
last ten years. During predesign studies for sewer 
outfalls, current meters were depioyed for periods of a year 
DEDO re near Santa Cruz [Brown and caldwell, Inc., 1978], 
„> Pajaro River [Environmental Research Ccnsultants, Inc., 
КОО and tha Salinas River [Znginsering-Science, Inc., 
1977]. In general, these studies suggest that «he net flow 
of water is northward along the coast, with some dependence 
en the iccal iiurnal wind. 

Tidal forsing has not been examined closely in any CL 


the aforementioned studies. However, an xplicit numerical 


с 


model was employed by Lazanoff (1971) to study the tidal 
circulation within the bay. Although field observations 
Meaweated that the primary driving force for circulation 
derived from »ceanic currents, the tide- and wind-driven 
model did predict correct sea-surífíac2= slevations and current 
phases and directions for the short time periods over which 
it coulá be run. Current magnitudes appeared too large near 


eaS” tal boundaries. 
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More recentiy Bretschneider (1980) analyzed the effects 
of various oceanographic conditions on the sea-level changes 
cbserved at Monterey. Variations in geostrophic current 
flow, atmospheric pressure, sea-surface temperature, and 
meridional wind stress were shown to correspond to observed 


variations in sea-surface elevation at Monterey. 
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II. DESCRIPTION OP IHE NUMERICAL MODEL 


— са: олы ли —— ა лын AGB Ee 


A. HYDRODYNAMIC THEORY 


The numerical model of Leendertsa (1967) relies on the 
basic equations that describe conservation of momentum and 
mass for incompressible fluid motion. In the Cartesian 
coordinate system of the model, with x- and y-axes embedded 
in a horizontal f-plane tangent at the origin to the undis- 
turbed sea surface (the datum) and with the z-axis oriented 


upward, these well-known =quations are: 


ბს ბს ბს ბს 185 


Se MS аа шшш, ш” >, (2.1) 
ra + ape + vet + wet _. (2-2) 
E + u; + თა = => ЗА (2.3) 

Rem = 0 (2.4) 


The variables u, v, and w are components of velocity 
parallel to the x-, y-, and z-axes respectively, p 15 pres- 
sure, and p is density. The applizd forces per unit mass 
(წ. ) represent effects of the Earth's rotation, the Earth's 
Gravitation, viscous and turbulent stresses in the fluid, 
and astronomical tides. 
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These equations are simplified by making assumptions 
appropriate to the examination of long-period forcing in a 
two-dimensional, shallow field. Detailed derivations may be 
found in Leendertse (1967) and Hart(1976). A more general 
development of shallow water equations may be found in 
Csanady (1982). The necessary asSumptions are summarized in 
the following paragraphs. 

First, because a coastal sea or estuary is generally 
shallow relative to the horizontal scale of motion, the 
vertical velocity 1s assumed small relative to the hori- 
zontal velocities. Therefore, both convective-inertia terms 
and rotational effects that involve the vertical velocity in 
equations 2.1-2.3 may be neglected. 

Second, the equations 2.1-2.4 are averaged to model 
fluid motions with periods greater than those of short- 
period turbulent motions. 

Third, the hydrostatic approximation is made by analysis 
of əquation 2.3. The vertical component of the rotational 
effect пау be considered negligible (of the order 10-2 
cm/s?) and so may vertical stress-gradient effacts (10-3 
cm/s? according *o Csanady, 1982). Furthermore, since mean 
vertical velocties are unlikely to be greater than 10 cm/s, 
over sufficiently long time periods (>103 s or 15 minutes) 
the total vertical acceleration will be of similarly small 
order. Neglecting for the moment tidal effects, an expr=s- 
Sion for pressure may ke obtained by integrating the 


remaining terms over depth: 


T) 


р = p t9 , 062 


(2.5) 


In this expression, n is the sea-surface elevation and р. 15 
the atmospheric pressure at the sea surface, both functions 
Of x and y. The gravitational acceleration, g, is assumed 
constant and equal to its mean valua at the undisturbed sea 
Surface. 
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Equation 2.5 permits computation of pressure gradients, 
S6e/ 6x and бр/бу, from horizontal gradients in sea-surface 
elevation and density. Gradients in atmospheric pressure 
may be neglected since their effect is small relative to the 
turbulent stress induced at the surface by the wind. Іп the 
numerical model, tidal forcing is accomplished by generating 
gradients of sea-surface elevation rather than by attempting 
tc simulate directly the astronomical forces tha- cause tke 
tides. 

Fourth, the BoussinesgJ approximation is made, in which 
the influence of vertical density variations is assumed -o 
be negligible. For this to be true tha arsa to be modelec 
must be vertically well-mixed, an assumption that is not 
generally applicable. Although the version of Leendertse's 
mcdel used in this study requires this assumption, some 
compensation for density stratification might be made by 
modifying the model to integrate estimated values for the 
vertical density variation over the depth at each point to 
obtain «he additional density-induced sea-surface elevation 
Mesenecy, 1982]. 

Fifth, the mean visccus-shear stresses of the fluid are 
assumed negligible, leaving only the turbulent stresses 
(Reynolds stresses) at boundaries within and external to the 
fluid to be formulated. Of these, Sharp density gradients 
within the fluid are neglected as a source oi stress. 

Cissed lateral boundaries are considered by applying the 
coastal boundary condition that velocity into the boundary 
is zero. Two other boundaries are considered, the sea 
surface and the bottcm, and algorithms for modeling stresses 
on these boundaries must prepared. 

Sixth, since interior stresses resulting from sharp 
density boundaries are assumed negligible, equations 2.1 and 
2.2 may be integrated vertically to provide implicit expres- 


sions for horizontal velocities averaged over depth. 
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Beplying the kinematic boundary condition at the free 
surface and at the bottom (assumed impermeable), equation 
2.4 may also be integrated verticalily using the Leibnitz 
mee ctf integration. Three equations implicit in three 
unknowns ars then available for the vertically integrated, 
horizontal velocities U and У апа the free, ses-surface 


elevation n: 


80 $0 OU. он В (2.6) 

Tag Т VE = 9 rd а — Fay) /P (h+n) 

SV SV м 3 

Da s VET = -а 2 £U + წეს ჯაი” (2:7) 
$n „ Sthmu , У _ე (2.8) 


6% 5х 5у 


In these equations f is the Coriolis parameter, h is the 
depth, Р 15 the surface friction stress due to wind, and S 
is the bottom friction stress. U and V are mean horizontal 
velocities over the water column. This sinplification 
results in a two-dimensional model with which patterns of 
circulation and sea-surface elevation may be quantitatively 
determined. 

Finally, since bo--om stress depends on the fluid 
velocity, a well-known quadratic model is assumed so that 
the stress term may be incorporated directly into 
Leendertse's computaticnal scheme. Tha formulation for 
bottom stress is: 
* 


F = ogU(U^«v^) ис? (2.9) 


Вх 
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= იძV(0?+V?) ус? u) 


Fay 
The empirical Cheézy coefficient, C, may be computed in any 
of various ways and must be specified for the area to be 


mcdeled. 


В. STRUCTURE AND COMPONENTS 


A derivation of the numerical model, a discussion of its 
computational stability, and a program listing are given by 
Leendertse (1967). Key features of the computational scheme 
and the computer program used during this study are 
presented hers. 


1. Computational Schene 


In the numerical model, equations 2.6-2.8 are 
approximated with a finite-differencing scheme extending 
over two time levels, each a half time step. At the first 
half time step, t+1/2, the velocity U(t+1/2) and the sea- 
surface elevation n(t+1/2) are computed implicitly and the 
ШЕ Осісү V(+t+1/2) is cemputed explicitly. t the second 
halz time step, t+1, the velocity V(t+1) and sea-surface 
elevation n(t*1) are computed implicitly and Ult+1) is 
computed explicitly. 

Computations are spatially controlled by a uniform 
grid of squares laid over the f-plane (Figure 2.1). Depths 
relative to the undisturbed sea surface, hare taken to be 
mean lower low water (MLLW) , must b= supplied for the 
Ccrners of each square, values of horizontal velocity are 
computed at the centers of the sides of each square, and 
values of sea-surface elevation are computed at the center 
of each Square. Wind-stress and bottom-friction factors 
must be specified cr computed at the centers of each square. 


This staggerei-grid is basic to the spatial realization of 
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Figure 2.1 Staggered Grid of the Wumerical Model. 


Leendertse's finite-di£ferencing scheme in that he mean 


MONO city into or out of each Square is used to compute the 
change in sea level within. The grid also permits the 
coastal constraint cf zero transport perpendicular to 


sea/land boundaries. 

ВЕ the first half tine step, implicit computa- 
tions proceed row by row from left to right and explicit 
computations proceed column by column from the bottom to the 
top of the grid. During the second half time step this 
process is reversed, "centering" the results in space as 


well as time. 
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The computer program used to model Monterey Bay is a 
modified version of a Fortran program developed by Hart 
(1976). It includes provisions for modeling wind stress, 
Steady flow at boundary points, and overflow at boundary 
pcints, none of which were available in Leendertse's orig- 
inal listing. To increase the flexibility of the program, 
additional modifications have been made during this study. 


These include: 


e Introduction of date and time computations 
to permit the program to search time-coded 
files for data items required at each half 


time step. 


* Direct computation of bottom-friction 
factors from an average bottom-type parameter 


or from a grid of bottom types. 


• Creation of an int2ractive subroutine to 
start the program by prompting the user for 


parameters critical to each run. 


e Output of time series of currents and 
Sea-surface elevation for up to nine points 


in the model grid. 


e Further modularization of the model's 


functional components. 


The core of the program is subroutine MODEL, which 
contains the finite-differencing scheme of Leendertse. 
Other subroutines serve auxiliary functions: Interactively 
starting the run, acquiring both constant and time-dependent 
daa”, Specifying conditions at boundary points, specifying 
numerical models for wind and bottom stress, and supplying 
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results in various output formats. An outline of the 


program is presented in Table I. 


TABLE I 
Components of the Numerical Model 


MAKEGRID Generate depth and computation-coatrol grids. 


BAYMODEL Main program to control each run of the model. 
START Interactively input run-time parameters. 
DIMENS Input constant data for grid of arsa. 


FIND Locate water sections to be modeled. . 
INVAL Initialize variables everywhere in grid. 
INCURR [nitialize currents as desired. 


MODEL Multioperational finite-difference scheme. 
CHEZY Supply Ghezy coeffrcrspt at each grid point. 
WIND И wind stress term at each grid point. 
RESULT Compute output values in desired units. 





© 
<j 
mj 
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OPEN Specify sea-surface elevarcions at open bounds. 
STEADY Specify currents at open bounds and rivers. 
Specify overflow currents at boundaries. | 
OVFLD Specify overflow threshholds. | 
INTIDE Set tide from a time-coded file. | 
INWIND Set wind from a time-coied file. | 
? | 
HEADS Jutput headers for each ӘШЕРШ- rb | 
PRINT Output results for desired times. | 
PLOT ВЕЧЕ ლ) == needed fer graphic display. | 
SERIES Output results for specific points. | 
PTGRID Utility — of input gridded data. | 
CALEND Utility for number of days in month. 
ADTIME Utility to increment time. 
BAYPLOT Provide graphic presentation of results 
| ELEVCOMP Compare elevation series with observed data. | 
| CURRCCMP Compare current series with observed data. | 
| 
| | 
ocon a Е ———————————tÀ 


Some subroutines, such as OPEN, must be prepared 
specifically for the arsa to wnich the model is applied, 
while others, such as MODEL, should not be altered. 


Modularization of the program permits the user to readily 
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change the appropriate functions when adap-ing the model to 


different coastal areas (and computers). 


C. USE AND ADAPTATION 


When applying the numerical moi2l to a specific area, 
certain requirements concerning input data and program modi- 


fication must be met. These are discussed below. 


1. Input Data Reguirements 





Each run of the numerical model requires specifica< 
tion of start and end times, time-step length, and an 
interval at which results must be output. For experimental 
(as opposed t» operational) use, other things may be speci- 
fied: Pcints at which series output is desired, the type of 
output, and onission of certain terms in the hydrodynamic 


equations. 

Input values that are uniqu2 to the area to be 
modeled and that do not change from one run to the next must 
also be supplied. These values are most conveniently stored 


in a separate file. They include: 
и Рос“ + ов title and central latitude. 
e Dimensions of the grid. 
e Depths for each grid corner. 


e Ccntrol numbers for each grid square 


(land=0, water=1, over flow=2). 


e A general bottom-friction parameter or a 
bottom-type indicator for the canter of each 


grid square. 


e Number of tide stations supplying data for 


boundary conditions. 


e Number of wind stations supplying data. 
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Selection of the size of grid to be used is limited 
by the size of the area to be modeled and the available 
virtual memory and CPU time of the computer. Since depths, 
ei ction and wind factors, output data, and two half-time- 
step values for each velocity component and the sea surface 
elevation must be available for each grid square at all 
times, at least 12 arrays must be dimensioned according tc 
the grid size and survey area. The maximum dimensions of 80 
by 80 used in this study required close to 1 megabyte of 
virtual computer memcry and approximately 0.5 s of CPU time 
per time step on the IBM 3033 mainframe computer. 





2. Subroutine Modifications 


Stresses at the bcttom and surface must be modeled 
in the subroutines CHEZY and WIND. Since Lesndertse's model 
already assumes a quadratic formulation for bottom friction, 
only the Chezy factor, C, need be provided by the subroutine 
CHEZY; however, many empirical techniques exist for 
computing the factor, most of which rely on a description of 
the bottom type. The user may select the technigue which 
best applies. 

Similarly, the user must program a wind-stress 
fermulaticn in the subroutine WIND. Values for wind speed 
and direction are obtained from time-coded data sets using 
the subroutine INWIND. 

Subroutines OPEN, STEADY, OVFLO, and OVFLD supply 
time-varying values for sea-surface elevation, currents, and 
overflow conditions at both open and closed boundary points 
in the grid. In OPEN, ап ailgoritha must be provided to 
compute the variation of sea-surface elevation along the 
open boundaries of the grid. The necessary tidal amplitudes 
are obtained from time-coded sets of data for established 
tide stations, using subroutine INTIDE. STEADY permits 


currents to be assigned to individual grid points, 
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Everridrng the computed currents. Ап initial current field 
may be enterel using subroutine INCURR. Finally, OVFLO and 
OVFLD permit conditions of flooding to be ascsrtained and 
modeled at grid points assigned tha control value 2. 

When implementing the model on various computers 
and for various purposes, modifications may be necessary in 
the input/output subroutines START, DIMENS, HEADS, PRINT, 
PLOT, and SERIES. RESULT may also De modified to compute 
additional quantities of interest, such as horizontal 


transports. 


3. Computer Implementation 





For this study, the numerical model was implemented 
on an IBM 3033 mainframe computer at the Naval Postgraduate 
School. Only a few minor language changes were required 
before the system's Fortran 4 compiler could be used on the 
program originally supplied by Hart. Subsequent modifica- 
tions were made and all jobs were run from remote terminals 
under the School's interactive +tima-Sharing systén. 

The system made possible the davelopment of several 
programs that facilitated preparation of data for input to 
me noel and production cf graphic output. Especially 
useful among these were: MAKEGRID, a program that generates 
Euout depth and computation-controel grids from digital 
bathymetric data already available for the area, simplifying 
the otherwise laborious task of creating a grid on chart 
overlays; programs that generate predicted tidal amplitudes 
from constituents, or from data supplied in the NOS Tide 
Tables; and, ELEVCOMP and CURRCOMP, programs that plot 
time-series output frcm the mcdel ayainst observed data from 
the same time period for verification of model accuracy. 
Although the programs themselves may not be transferable to 
other computers, supplying similar auxiliary software 


together with the model (or even incorporating the 
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algorithms into the model program) greatly enhances the 
ready application of the model to other coastal areas. 

Data necessary to running the numerical model were 
stored in computer files distinguished by type. All 
constant, gridded data were stored in one file while time- 
varying data were stored in separate files by type and year 
(for example, MONTEREY TIDE76). In this way, a new file cf 
input data did not have to be created for =ach run of the 
model. The sources and selection of input data are 


discussed in sections 3.B and 3.C. 
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III. APPLICATION OF THE MODEL TO MONTEREY BAY 


А. VALIDITY OF ASSUMPTIONS 


The applicability of the numerical modal to Monterey Bay 
was checked by examining the assumptions outlined in section 
2.А. The assumption of nagligibl2 vertical velocity was 
confirmed for tidal forcing by noting that the maximum depth 
of the area modeled is much léss than the wavelength of the 
semidiurnal tide (3 km << 7600 km). Other, horizontal, 
forcing conditions of currents and wind were applied as 
steady-state phenomena in the model. In addition, since 
this study concerns large-scale fluid motions over “іле 
periods of 15 minutes or more, shorct-period turbulent 
effects and vartical accelerations were neglected: the 
hydrostatic approximation holds. 

Although Monterey Bay 1s not vertically weil-nixed 
(Scott, 1973], in depths of a few hundred meters or isss the 
difference in dynamic height between that of the assumed, 
homogeneous density profile and that of a nore typical 
Exe 1s less than 1 cm, which is negligible for the 
purpose of this study. In depths of 1000 m or more, the 
effect is mor? significant (several centimeters); however, 
Since such depths occur outsid= tha bay proper, the effect 
of density stratification was ignorad and horizontal veloci- 
ties were averaged over depth to obtain a general picture of 


eırculation in the bay. 


B. CONSTANT INPUT 


The numerical model was applied to three differen- grids 
covering Monterey Bay (Figure 1.1 and Figures 3.1 - 3.3). 
Grid A, a small-scale, 1- or 2-km grid, 80 by 80 km, was 
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designed to permit the introduction of offshore, non-tidal 
currents as a steady forcing condition and to place the bay 
far enough away from the three open boundaries 50 remove 
their associated spurious sffects. Since this grid was 
particularly vulnerable t> numerical instabilities in the 
model, two other grids were devisel. Grid B, a large-scale, 
1-km grid, 23 by 50 km, covered the bay and reduced the 
number of open boundaries to one. Srid C, a 1-km grid, 40 
by 72 kn, covered both the bay and sufficient area to permit 
offshore, non-tidal currents to be introduced. All griäs 
were skewed 209 east cf north to align the boundaries 
perpendicular to the tidal forcing conditions. The dinen- 
sional and constant data incorporat2d into these grids are 


discussed below. 


1. Time and Space Dimensions 


A time step of one hour was chosen to permit assess- 
ment of the model over pericds of several days without 

ating extensive us2 of CPU time. Normally the model 
Exculd run for 12-24 hours (one tidal cycle or nore 
depending on the tidal phase differences between varicus 
Bes or the area| to establish reailistic conditions of 
current and sea-surface elevation throughout thes area (Hart, 
EN]. The one-hour interval provided a sufficient number 
СІ дата values fcr ccmpatison with hourly or half-hcurly 
observaticns Jf sea-surface elevation and currents. 

Since computed values are offset in the staggered- 
grid scheme of the mcdel, use of a relatively small grid 
Siz2 in regiens of steep bathymetric relief is important to 
model accuracy. With the constraint on array dimensions in 
Mind, the smallest grid size possible (1 km?) was generally 
selected tc permit the greatest spatial resolution for the 


area of concern. 
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Depth Contours for Grid A. 


Figure 3.1 
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Depth Contours for Grid B. 
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Depth Contours for Grid C. 


Figure 3.3 
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Digital bathymetry was provided by th= NOAA National 
Geophysical Data Center, where depth data from past hydro- 
graphic surveys are archived for most coastal areas of the 
United States. The depths were positioned by latitude and 
longitude in à 36-sec grid and were referenced to a mean- 


lower-low-water (MLLW) datum. 
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Figure 3.4 Locations of Bathymetric Data in Grid A. 


To conputer-generate a grid of depths for the model, 
tha program MAKEGRID was prepared. The program first 
projects the bathymetric data onto the grid coordinate 
system (for example, grid A), which is a modified-transverse 


Mercator projection skewed about a specified origin 
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(Figure 3.4). Depths are then interpolated at the corners 
of each grid square. 

code = аге raferenced to MLLW and a straight- 
forward correction to mean sea level is not possible, the 
datum for sea-surface elevations computed by the model was 
taken to be MLLW. All input tidal amplitudes ware likewise 
referenced to MLLW. 

omenez i dced depths a computation=cóntrol grid 
was automatically generated by assigning 1's to all water 
squares and 015 to all land squares (assigned dummy eleva- 
tions in the lepth grid). Both grids could be altered if 


necessary before their use in the model. 


To model bottom stress, the empirical Manning equa- 


Eon for the Chezy coefficlent, C, was used: 


А м (3:91) 


The coefficient is а function of depth, h, instantaneous 
sea-surface elevation,n , and the Manning factor, M, which 
describes bottom roughness. M increases with bottom rough- 
ness. Clean and straight natural river channels typically 
require M=0.025 to 0.030 m/s, while winding chanrels may 
require M=0.033 to as high a value as 9.15 m/s in Very 
weedy, overgrown areas (vo. 99, Henderson, 19661. 

Although bottcm stresses may be modeled as a func- 
tion cf the bottom type or texture in each grid square (thus 
шј па input of a bot<om-type grid), this option ‘was not 
exercised for Monterey Bay. Over the large area covered, 
the variation in depth from square to square is likely to 
have more influence than the relatively small variations in 
bortom roughness that occur in Monterey Bay. Following 


Spaulding and Beauchamp's study of a coastal sea (1983) and 
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after some experimentation (section IV.A), a constant value 


cf M=0.04 m/s was used throughout the bay. 


C. TIME-VARYING INPUT 


The forcing conditions of tide and wind were accessed 
from time-coded files. Tidal amplitudes were appiied only 
along the open boundaries of the model, whereas wind stress 
was applied cver the entire grid. The sources and applica- 


tion cf these data are discussed in this section. 


At the start of a run of tha model, sea-surface 
elevations were approximated by assigning a coastal tidal 
amplitude at the starting time to every point in the grid. 
Tha tidal amplitude at Monterey was used for this purpose. 
This approximation is Suitable for Monterey Bay since the 
narrew continental shelf and the absence of any barrier 
isiands permit the tidal wave to propagate relatively 
rapidly through the area. 

À zero velocity was initially assigned to each grid 
Beene, except for runs including an offshore, non-tidal 
current; in these cases, the steady-state current velocity 


was initially assigned to offshore grid points. 


2. Boundary Tidal Amplitudes 


A major factor in the successful application of he 
nunerical model ws the provision of suitable tidal forcing 
conditions along the open boundaries. Tidal amplitudes are 
predicted in the NOS Tide Tables 1976 for four stations near 
and within Monterey Bay: Ano Nuevo, Santa Cruz, Monterey, 
and Carmel. Tidal-constituent amplitudes and phases are 
available for Monterey and Moss Landing (Appendix A). Since 


1963, continuous observations of water level have been made 
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at Monterey. A two-year series of nearly continuous obser- 
vations was made at Moss Landing from 1976 through 1977. To 
obtain amplitudes along the boundaries of the model grid, 
coastal values such as these must b2 extrapolated. 

Several attempts have been nade to formulate the 
effects of continental shelves on the open-ccean tide 
(Clarke and Battisti, 1980; Gill and Porter, 1980; Munk, et 
al., 1970]. Because of the narrow continental shelf and 
bisecting canyon, Monterey Bay does not satisfy the assump- 
tions necessary to apply these formulations. However, to 
gain some insight into the effect of the extreme depth 
difference between the tide station at Monterey and the 
offshore boundary points of the model, a comparison was made 
between tidal constituents obtained at Monterey and at a 
pressure gage located in 3903 m of water offshore 
(Cartwright, ət al., 1979]. The results are presented in 
Table II. 

The coastal and pelagic phases clearly Go not corre- 
spond since the pelagic gage was located at some distance 
from Monterey (see Figures 1.1), but *he agreement of the 
amplitudes suggests chat the aforementioned depth difference 
has little effect. Inthe absence of any more certain 
method for extrapolating tidal amplitudes, the values at “he 
coastal station were applied directly along a line of 
constant phase extending out from shore. 

Examination of cotidal/cophas? charts by Munk, et 
БИЕ (1970), Luther and Wunsch (1975), and Parke and 
Henderschott (1980), suggests that, in the vicinity of 
Monterey Bay, the tidal wave propagates nearly parallel to 
the coast. The model grids were, therefore, skewed in such 
а мау that the open boundaries wer2 perpendicular ог 
parallel to the coast. At each time step the forcing ampli- 
tude was made to vary directly with the tide at Monterey all 
along the southern boundary, with the tide at Ano Nuevo ali 
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T ABLZS II 
Comparison of Pelagic and Coastal Tidal Constituents 
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along the northern boundary, and with linearly interpolated 
values between the two along the western Doundary. 

The tidal values used to sat the boundary conditions 
may be interpolated from the NOS Iiis Tables 1976, computed 
meom Constituents [Schureman, 1950], or taken directly fron 
observed data. The last was preferred since the first two 
predictive techniques cannot take into account atmospheri- 
cally forced or anomalcus changes in sea level, such as 
storm surge. Howsver, as mentioned previously, observations 
were available only at Monterey ani Moss Landing. Some 
experimentation was necessary to model the phase lags 


between Monterey, Santa Cruz, and Ano Nuevo (section IV.A). 
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3. Boundary currents 


Two types of flow may be imposed at the boundaries 
of the numerical aodel. First, the discharge of rivers 
along otherwise closed boundaries may be represented as a 
vertically averaged current velocity assigned to the appro- 
priate coastal grid point during each time step. The mean 
annual discharge of all major streams and rivers entering 
Monterey Bay is 1.85 x 109 m3/day [Broenkow and Smethie, 
1978 ], which amounts to a vertically averaged current 
velocity of only 0.2 cm/s were all rivers to enter at one 
point. As a result, the river inflow was judged negligible 
БӘР this application. 

A second type of flow, currents due to non-tidal 
effects, may be imposed in offshore regions of the model. 
The narrowness of tne continental shelf near Monterey Bay 
leaves the bay particularly open to forcing by large-scale 
oceanic currents. Previous studies of the area suggest that 
Бол currents are an important force driving the circulation 
DE the bay [Lazanoft, 1971; Garcia, 1971; Bretschneider, 
1982]. The presence of an offshore current was simulated by 
HESENUMeng initial velocities to the offshore portion of 
grids A and C for some runs of the model. The convectivs- 
inertia *erms in che numerical model propagate the current 
ши Шшепсе 1п:5 the inshore portions of the grid. A north- 
ward current 9f 25 cm/s (0.5 knots) was assigned. This 
va_ue was proposed by Scott (1973) as a simple, steady-state 


model for the offshore circulation. 


и. Wi 


nd 


Wind stress, F., was parameterized within the numer- 
L 


ical model by the familiar quadratic law [Dronkers, 1964]: 
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Pi = (2.6x10°")p_ w|w|/p(h+n) (3.2) 


W represents the wind velocity vector, р. is atmosphetic 
pressure, p is the mean density of the water, h is the 
depth, and n is the time-varying saa-surface elevation. The 
model permits input of wind direction and amplitude as a 
ЕС 109 condition over the whole field of the grid for a 
specified ranje of time steps. Monthly distributions of 
wind at Santa Cruz and Moss Landing for the period May, 
1976, through May, 1977, were obtained from the Santa Cruz 
Wastewater Facilities Planning Study [Brown and Caldweil, 
Inc., 1978]. Average and maximum values for the wind were 


applied to sone runs of the model (ses section IV.B). 


D. DATA FOR COMPARISON 


The numerical model was calibrated by comparing modeled 
sea-surface elevations and current velocities at specific 
grid points with observed values at the same locations. The 
comparison process was limited by a paucity of suitable, 
long-term water-level and current-m2ter observations for 
Monterey Bay. Water-levei data are available only fcr 
National Ocean Service tide stations at Monterey and Moss 
Landing. The primary sources for current-meter data are 
predesign studies conducted for the emplacement of sewags 
outfalls near Santa Cruz, the Pajaro River, and the Salinas 
River, but only data for Santa Cruz and the Salinas River 
could be obtained. In some fortuitous instances, both 
water-level and current-meter data were collected concur- 
rently (Figure 3.5). The period July-August 1976 yielded 
sufficient data for comparisons at the two water-level and 
at three current-neter stations; their locations are plottéd 
on each model grid (Figures 3.1 - 3.3). 
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Figure 3.5 Observíng Periods for Comparison Data. 


1. Mater-1 


yo 
t1 


vel Observations 


= 


ater-level observations have been made nearly 
continuously since 1963 at NOS Tide Station 941-3459 on the 
seaward end of Municipal Wharf 2 in Monterey. The float- 
type tide gage is located in 5.2 m of water. Recorded times 
are accurate to within 6 minutes and heights are resolved to 
Вся | Bretschneider, 1980]. Digitized hourly heights for 
the period 11/8/73 to 3/2/83 were obtained from the National 
Ocean Service, Tilal Datums Section N/OMS123. The heights 
were correctei to MLLW and converted from feet to meters. 

In addition to providing comparison data, these observations 
were used to determine boundary amplitudes for some runs of 
the model. 

At Moss Landing water-level observations were made 
for 20 months as part of the California Marine Boundary 
Program [National Ocsan Survey, 1981]. NOS Tide Station 
941-3616 was a float-type tide gage located at the seaward 


end of the Moss Landing Ocean Pier in 9.1 m of water. 
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Digitized hourly heights were obtained for the entire period 
5/9/76 to 1/10/78 and processed as for the station at 


Monterey. 
сос теле 


As part of the Santa Cruz Wastawaters Facilities 
Planning Study, a current-meter station was located 1 mile 
EBrcshere of Terrace Pcint in 30 m of water { Brown and 
Caldwell, Inc., 1978]. Two AMF Vector Averaging meters were 
installed at 9- and 15-m depths for the periods June to 
November, 1976, and January to May, 1977. The only data 
that could be obtained for comparison purposes covered July 
and August, 1976, az the 15-m depta. The data included 
7.5-ninute averages cf current speed and direction, 
exoressed as a pair cf orthogonal velocity vectors. 

Two current-metar stations wer2 occupied approxi- 
Mazely 1 nautical mile north aad south of the Salinas River 
đuzing oceanographic investigations for the Monterey 
БЕЛЕГ 501а Маса Pollution Control Agency 
NESUSDHe2sring-Scienca, Inc., 1977]. at 2ach station, two 
ducted-impeller-tyoe meters were installed at 9 and 15 п for 
the overall pericd January, 1976, tə January, 1977. Current 
speeds and direction were averaged at 30-minute intervals 


and expressed as a pair of orthogonal velocity vectors. 
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IV. RESULTS 


A. COMPARISON OF MODEL RESULTS WITH OBSERVATIONS 


Comparison of the model with observations was used both 
to "fine tune," or calibrate, the numerical model, and to 
assess its general validity. The 3ffects of varying input 
constants such as the size and resolution of the grid, the 
time step, and the bottom-friction coefficients were consid- 
er2d. In additicn, schemes for determining boundary tidal 
amplitudes and for including non-tidal current fields were 
tested in an effort іс match observed elevaticns and 
cuzrents as closely as possible. 

Application of the model to grii A revealed apparent 
numerical instabilities that caused overflcw in the computa- 
tions after as few as 1.25 days (39 time steps). The sudden 
oscillations in sea-surface elevation at Monterey were due 
to propagation into the bay of extreme amplitudes and 
тепе generated in the offshore portion of the grid 
(Figure 4.1). The overflow condition was unaffected by 
вазе the resoluticn of the griji from 2 to 1 km, but was 
very sensitive to changes in the phasing or the tidal ampli- 
tudes along the open boundaries. Under the premise that tne 
presence of three open boundaries 2nhanced instabilities, 
grid B {with one open boundary) and grid C (with two open 
boundaries) were subsequently used iuring the comparison 


process. 


1. Sea-surface Elevation Comparisons 





A run of the model for twelve days at the one-hour 
tine step produced scme agreement between nodeled and 


observed sea-sur face elevations at Monterey and very good 
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Figure 4.1 Sea-surface Elevations at Monterey, Grid A. 


agreement at Moss Landing (Figures 4.2 and 4.3). Varying 
Ens Manning bottom-friction factor improved this result 


(Table 3.1 and Figures 4.4 and 4.5). 


ыы... д 


TABLE III 
Effect of Various Manning Factors 





| | 
| 

| M (m/s) = • 03 • 04 2915 26 . 10 | 
| Monterey . 7.8 6.8 6.3 5.8 4.8 | 
| Moss Landing 4.6 15 4.4 4.3 4.0 | 
Table values are the RMS errors in centimeters | 
for ccmparisons between modeled and observed | 
sea-surtace elevations. | 

1 | 
–--““““-__'-____-“-“_”“”“”რ”ნწრნჩნ eK mM MEME mE с UEM amp ———- ЕСЕЈ 
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The seemingly unrealistic Manning factor of 0.10 ıYs 
may serve as 2 description of the bottom roughness over the 
large area of each qrid square, if roughness is thought of 
in terms of the steep slope that is otherwise not accounted 
for in the model. The higher Manning factors did not, 
however, improve the results of current comparisons; they 
served primarily to damp out noise. 
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Figure 4.2 Elevations at Monterey, ^tz3600 s. 
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An oscillation that appeared forced by the time step 
was evident in the modeled curves. It was especially 
evident at Monterey and when a shorter, 15-minute time step 
was used (Figures 4.6 and 4.7). 

The noise may be the result of applying observed 
tidal amplitudes as the boundary forcing condition. The 
cbserved water levels, digitized hourly, пау include jumps 


and/or contaminating frequencies due to the recording 
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Figure 4.3 Elevations at Moss Landing, Atz3600 s. 


instruments. Linear interpolations required between rhe 
hcurly amplitudes for each half time step may have exagger- 
ated the instrumental effects. To better judge the use of a 
Shorter time step, raw water-level observations, usually 
made at a 6-minu<-se interval, shouid be applied to the model. 
Spectral analysis of the modeled and observed 
curves, in addition +o reflecting their general agreement, 
reveals spurious frequencies generated by the model at 
Monterey (Figures 4.8 and 4.9). Tha incoherent frequencies, 
meen are als> found in the spectra for currents at Santa 
ezuz (Figure 4.14), correspond to apparent periods of 3.0 
=>) 2.2 hours. These periods are longer than the 1-hour, 


fundamental seichne pericd for Monterey Bay (Lynch, 19701. 
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Figure 4.7 Elevations at Moss Landing, At=900 s. 
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Figure 4.8 Spectra, Elevations (SE) at Monterey. 
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Figure 4.9 Spectra, Elevations (SE) at Moss Landing. 
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The phenomencn of time-step-linked oscillations has 
been experienced by other investigators (Chiang and Lee, 
1982], who found that generating a nathematically smooth 
function from the observed data provided more suitable 
amplitudes for forcing the model. A smooth amplitude func- 
tion was obtained in this study by summing the tidal 
constituents for Monterey. Applying boundary conditions 
based on these predicted tides gradually reduced the noise 
in the model (Figure 4.10). Not surprisingly, however, the 
predicted elevations did not provide good agreement with the 


observed water-leveis. 
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Figure 4.10 Elevations Using Predicted Boundary Amplitudes. 
The sensitivity of the model +o open-boundary condi- 


tions, already noted in the case of grid A, was demonstrated 


by compazing results for different phase lags along the one 
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boundary of grid B (table IV). A constant amplitude along 
the boundary (no phase lag) was found to be most successful, 
though the NOS Tide Tables 1976 predict that Monterey lags 
Santa Cruz by 6 minutes. 


( 
| 
| 
| 
| 
| 
| 
| 
| 


TABLE IV 
Effect of Boundary Amplitude Phasing on Grid B 


Phasing 6+ 0 6- 
Monterey 7.9 7.4 7.8 
Moss Landing 5:53 4.6 5:0 


Table values are the RMS errors in centimeters. 
Phasing is Monterey minus Santa Cruz, in minutes. 


жымы as ————————— 


" 


2% Current Comparisons 


Тһе modei-generated current values agreed poorly 
with the 15-m depth observations at all three locations 
where comparisons were made (Figures 4.11-4.13). In making 
the comparisons, current vectors from the model output and 
from the cbserved records were resolved into eastwazd ard 
ncrthward components, taking into account their respective 
skewed coordinate systems. The mcdeled current components, 
particularly near the Salinas River nouth, were weak or 
non-existent. Clearlv, forces in addition *o tides were at 
work in generating the observed currents. 

For the currents near Santa Cruz the apparent rough 
cctrespondence in frequency was confirmed somewhat by a 
spectral analysis (Figure 4.14). The gradual increase іп 
the northward component of the modeled current may reflect 


leng-period variations in the current field. 
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In an effort to obtain better current agreement and 
to see the effects of currents associated with the conti- 
néntal slope, a steady offshore current was applied to grid 
С. The results near Santa Cruz were poor; inherent osciil 
tions in the nodel were amplified and agreement was not 
improved. The model's sensitivity tə open-boundary condi- 
eons Made this a difficult subject to pursue within the 
scope of this study. 


В. MODELED CIRCULATION OF MONTEREY BAY 


Although the tidally forced numerical mod2l did not 
reproduce obsarved currents a+ the comparison locations, it 
should nevertheless provides an estimate of the barotropic 
tidal circulation of Monterey Bay. A general view of circu- 
latory patterns may be obtained by sxaminiag the modeled 
sea-surface elevation and current fields for a 24-hour 


period. 
ee таат у Forced Circulation 


A tidally forced ¢levation-field series is presented 
in Appendix B. A small cscillatory structure in the 
southern bight ís consistent with the greater amplitude of 
noise observed a+ Monterey. In other respects, the eleva- 
tion fields generally show a clear progression of the tidal 
wave into the bay. 

A series of current-field plots, including th= volu- 
Metric transport associated with each vector, is presented 
in Appendix С. A total volume of 2.0x109 m3 appears +o be 
pumped across the boundary during the diurnai tidal cycle 
This barctropic result may be contrasted with the 109 m3 
pumped by internal +ídal mixing reported by Broenkow and 
Smethie (1978). 
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Figure 4.11 Currents near Santa Cruz. 
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Currents North of the Salinas River. 
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The current plots show generally weak (<5 cm/s) 
tidal flow into and cut of the bay, corresponding to periods 
ӨК + 1C6090 зпа ер. In theo southern part of the bay, a strong 
east-west current (up to 30 cm/s) appears just north of the 
Monterey peninsula. This jet is consistent with the strong 
currents experienced by divers in the area.! In th northern 
Bact Of the bay, a brcader current (10 cm/s) flows along the 
depth contours. 

Of especial interest is a Currant pattern that 
develops between Ano Nuevo and Santa Cruz on current plets 
made using grid A (Figure 4.15). The gyre, which the model 
predicts to have speeds ranging frou 2 to 10 ca/s, is 
consistent with the observations of Carter and Kazmierczak 
(1968) who noted a closed circulation in the area with 


similar speeds. 


2. Тійа11у Forced Circulation with Wind 








When an average wind of 3 m/s (7 mph) from the 
west-northwest was applied over th2 entire field of grid В, 
Am m 


the tidally forced circulation was unchanged. 


ct 
ის 
i а 
(Л 
O 


wand of 10 m/s (30 mph) from the west-nortiwesi 
produced essentially unchanged circulatory patterns. А 
12-day series cf modeled elevations at Monterey and Moss 
Landing under the same maximum wind yielded resultes nearly 


Шізісіса!і го those without wind. 


lPer conversation with Dr. E. C. Наде 


aderlie, Dept. of 
Oceanography, Naval Postgraduate Scnooi, 9/27/83. 
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V. CONCLUSIONS 


А. VALIDITY OF THE MODEL 


The numerical model has been shown by comparison with 
observations to provide reasonably accurate ssa-surface 
elevations. Although the modeled current doss not account 
for the total observed current at comparison stations, it 
probably accurately reflects the contribution due to the 
barotropic tile. Much of the remaiaing observed current may 
represent response to diurnal wind stress, offshore пог- 
mead. Currents, and/or forcing due to internal waves, “he 
last requiring a more complex, thrz=-dimensional model for 
further investigation. Application of observed winds ona 
timestep-by-timestep basis might bs a Zruitful avenue for 
further investigation. 

That numerical instabilities exist in the model has been 
hood by various authors [ Moe, et al., 1978; Benque, et al., 
1982], who have propcsed some inprz5vemen-s in the model's 
formulation. The flexibility and accuracy of the model 
Might be improved by further investigation of their 


suggestions. 


В. HYDROGRAPHIC SURVEY APPLICATIONS 


A numerical model such as that implement2d by this study 
can improve the process of correcting depths for changes in 
sea-surface elevation during a hydrographic survey. 
Advantages of the two-dimensional model cver simple extrapo- 
lative techniques or more complex, three~dimensional mcdels 
result from the model's relative sinplicity, flexibility, 
ability to operate in a real-time data collection systen, 
and ability t> compute sea~-surface 2lavations with suffi- 
Sent accuracy. 
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The tested model is a relatively simple FORTRAN program 
to implement and use, particularly with the interactive 
modifications made during this study. It could be further 
improved in this respect by adding an interrupt/restart 
routine to permit changes of constants, such as time step or 
outout frequency, during the course of a single computer 
run. The model is flexible, since it can be readily applied 
to various coastai areas by means of the gridding software 
developed during this study. 

To implement the model on a microprocesscr during data 
collection in the field, requires the availability of suffi- 
cient virtual storage capacity and CPU time to permit unim- 
peded computations. The virtual storage required depends 
upon the dimensions cf the grid; a nor?» economical use of 
arrays in the model program can reduce the requirement. In 
the real-time mode of operation, conputations should immedi- 
ataly follow boundary-amplitude updates at each half time 

сар to make afficient use of CPU «ime. At the conclusion 
of seach full tim? step the resulting, updated sea-surface 
elevations ar» then promptly available. 

The time step used depends upon the interval at which 
water-level observations are available from one or more 
locations suitable for establishing boundary amplitudes. Іп 
the real-time mode, presumably such data could be teleme- 
tered to the survey vessel at “he standard tide gage 
frequency of 6 minutes, permitting 2 model time step of 12 
Minutes. Sinc2 updated elevations can only be available at 
the conclusion of a time step (Figure 5.1), a $-minute lag 
exists “hat may be removed only by post-survey processing. 

Ànother factor in real-time operation cf tne model is 
the start-up time required. Тһе model should be calibrated 
me establish the validity of friction models and boundary- 
amplitude algorithms, preferably by comparing the output for 
several tidal cycles with observed tides at a location in 
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Figure 5.1 Time Step Lag During Real-time Model. 


the interior of the survey area. If historical data are not 
available, this could require several days of observation 
and analysis prior to the survey. 

The accuracy oí the model in coaputing see-surface 
elevations from tidal forcing alone has been estimated 
па this study as 4-8 cm (1 RMS 2rr9r). Survey require- 
ШЕШСЕ аге 30 < 9.14 * 0.005h cn, where o is the standard 
error ard h is the depth in centimeters [Mcblsy, 1982]. The 
depths at the tide gages in this study were about 8 п, 
permitting a 30 equal to 13.1 cm. This value was attained 
at Moss Landing (3 RMS error = 12 са) and, if “һе trouble- 
some noise could be removed by post-survey processing of 
water-level observaticns +o obtain a smooth tidal fercing 
function, it лау be generally attainable. 

The model itself requires further development, both in 
the application to Mcnterey Bay and in general. Further 
testing of the application to Monterey Вау should include 
the introduction of time-varying wind and oceanic currents, 
as well as forcing at a time step using tidal amplitudes 
observed at a shorter interval than the t-hour interval used 


in this study. Тһе model should also be applied to and 
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quantitatively tested against other coastal configurations 
and bathymetries te ensure that sufficient accuracy can be 


consistently achieved. 
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